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Abstract. We describe a cryogenic cavity-optomechanical system that combines 
Si3N 4 membranes with a mechanically-rigid Fabry-Perot cavity. The extremely high 
quality-factor frequency products of the membranes allow us to cool a MHz mechanical 
mode to a phonon occupation of n < 10, starting at a bath temperature of 5 kelvin. 
We show that even at cold temperatures thermally-occupied mechanical modes of the 
cavity elements can be a limitation, and we discuss methods to reduce these effects 
sufficiently to achieve ground state cooling. This promising new platform should have 
versatile uses for hybrid devices and searches for radiation pressure shot noise. 
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1. Introduction 

For decades the mechanical effects of light have been used to coax gas-phase atoms 
towards their quantum mechanical ground state of motion. Recently, experimenters 
in the field of cavity optomechanics have learned how to extend the mechanical 
effects of light to more massive mesoscale objects. A variety of microfabricated 
devices, with integrated optical and mechanical or electrical and mechanical resonators, 
have now been employed to backaction cool mechanical objects to near their ground 
state: Successful devices use lithographic electrical circuits [UEIIS], nanoscale optical- 
mechanical silicon resonators [4], or whispering gallery mode resonators [5]. These 
experiments combine cryogenic cooling of the mechanical object with backaction cooling 
of a specific mechanical mode by harnessing radiation pressure within a cavity. 

Parallel attempts have been made to cool mirrors or other dielectric objects to their 
ground state, and detect motion at quantum limits, within a canonical high-finesse 
Fabry-Perot cavity 0, [TJ El El QUI EU EEj. These devices are generally characterized 
by lower frequency, higher quality factor mechanics, larger-mass, and power handling 
capabilites closer to macroscopic cavities. Low-frequency mechanics combined with 
high quality factors enable longer absolute coherence times, and hence prospects for 
long-lived quantum memory elements. A larger motional mass is important for searches 
for gravitationally induced quantum collapse [13]. Another frontier in optomechanical 
systems is the realization of a quantum-limited continuous position measurement of a 
macroscopic object [TH UHl HH [T71 [181 HH ED]- This interest is in part motivated by 
large interferometers used for gravitational wave searches that will soon be limited by 
so-called radiation pressure shot noise (HI [15]. Thus far, operating optomechanical 
systems at high power has been one limitation to observation of radiation pressure 
shot noise, and Fabry-Perot cavities comprised of super mirrors and low-optical-loss 
mechanics are an excellent candidate for achieving the required intensities. However, 
despite the promising combination of optical and mechanical properties Fabry-Perot 
systems afford, solid-state mechanical resonators in Fabry-Perot cavities have not yet 
entered the quantum regime. One reason is that non-integrated cavities are difficult 
to make cryogenically compatible while maintaining optical alignment and vibrational 
stability. 

Fabry-Perot cavities are also particularly amenable to cooling high-tension, planar 
membrane mechanical objects [ID], which offer a promising route for creating hybrid 
systems for quantum information. In particular, there is considerable interest in 
combining optomechanical and electromechanical systems to realize mechanically 
mediated quantum state transfer between microwave and optical photons [2H E21 ESI 
EH llj- One of the successful platforms for electromechanics is a 10 MHz membrane 
coupled to a LC resonator, and hence one possible electro-optical transducer is a 
system where electrical and optical cavities are parametrically coupled to the same 
mechanical membrane resonator. However, no membrane mechanical system, and in 
fact no mechanical object with a frequency below ~ 70 MHz, has been brought to the 
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Figure 1. (a) Membrane-at-the-end cavity and simultaneous image of the intracavity 
optical mode and the SiaN4 membrane. The nodes of the (4,4) mechanical mode are 
indicated in blue, (b) Device design. Along the axial direction of the cavity from left to 
right are: Piezoelectric actuator (yellow), flat cavity mirror (transparent), membrane 
chip, silicon holder (green) with piezoelectric actuator (yellow), and curved cavity 
mirror (not visible). 

quantum regime with an optical platform [U [5] . 

Here we describe the details of a robust three-component Fabry-Perot cavity at 
cryogenic temperatures that incorporates silicon nitride membrane microresonators [10J. 
Our simple and near-monolithic design rigidly attaches the cavity mirrors and the 
membrane to a common base while maintaining the stringent alignment requirements 
of a high- finesse cavity [25]. With cryogenic pre-cooling, we can harness the extremely 
high frequency quality-factor products seen in higher-order modes of S13N4 membranes 
[26| [271 128| [21] , allowing, in-principle, ground state optomechanical cooling. By laser 
cooling this device in a 4 He flow cryostat at 5 kelvin, we achieve cooling of a few MHz 
mechanical mode to the lowest occupations yet achieved for an oscillator with this low 
a frequency, n < 10. The limitation to cooling is other background mechanical modes 
of the cavity structure; we discuss how to understand the effects of these modes based 
upon the transmitted intensity spectrum and future routes to eliminating these modes. 

2. Device and experiments 

Our cavity consists of two superpolished fused silica substrates with a high-reflectivity 
coating. The 100 ppm throughput mirrors can have down to a few ppm scattering- 
absorption, creating a cavity with a finesse up to 31,000. One mirror has a radius 
of curvature of 5 cm, and the other mirror is flat. For the work presented here, the 
mirrors are placed L — 5.1 mm apart and the membrane is placed 0.9 mm from the 
flat mirror. This creates a "membrane-at-the-end" cavity as shown in Fig. [T^a). An 
invar spacer connects the mirrors, the membrane, and the piezoelectric elements that 
translate the membrane and one mirror along the cavity axis (Fig. [TJV)). In contrast to 
early devices [25], the cavity design presented here allows direct access to the membrane; 
by pulling the entire central metal section vertically out of the cavity, the membrane 
can be switched out without disturbing the high-finesse cavity mirrors. 
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The cavity is assembled by aligning the elements at room temperature using an 
optical signal and then epoxying each piece in place in turn. We align both cavity 
mirrors to retroreflect a single fixed laser beam. The membrane is then inserted into 
the cavity, held temporarily on a 5-axis micrometer stage. The transverse alignment of 
the membrane is assessed by in-situ imaging of both the optical mode spot, formed by 
a 1064 nm laser source, and the membrane, illuminated with an LED source at 940 nm 
where the mirror reflectivity is low (FigJT](a)). Any tilt in the membrane plane relative 
to the optical axis leads to a displacement and distortion of the optical mode spot, 
which also varies with the location of the membrane along the optical standing wave. 
The membrane tilt is adjusted so the location and shape of the optical mode spot are 
not perturbed as the membrane travels over an optical wavelength. 

The cavity must be constructed in such a way to maintain alignment upon cooling 
the device down to cryogenic temperatures. To this end, we employ mainly low thermal 
expansion materials such as invar, fused silica, and the popular cryogenic epoxy Stycast 
2850FT. Further, a key element of the construction is symmetric, uniform thickness 
epoxy joints that do not result in relative tilt of the elements due to differential 
coefficients of thermal expansion, nor in excess localized mechanical stress build-up. 

The membrane inserted into the cavity is a square d = 500 /im Si3N 4 membrane 
from Norcada Inc. We use t = 40 nm thick membranes, as measured directly via 
ellipsometry. The membrane is suspended on a 5 mm square silicon chip with a 
thickness of 500 /xm. We focus on cooling either the (m,n)=(2,2) or (4,4) drum modes 
(where m, n are the number of antinodes in x, y) of the membrane. These modes have 
resonant frequencies of 1.6 MHz and 3.2 MHz and square mode sizes of 250 /im and 125 
/im respectively. These sizes can be compared to the cavity mode, which is typically 
measured to have an intensity profile with a 1/e 2 diameter of 180 /im at the membrane 
position. The mechanical mode frequencies are sufficiently large to theoretically cool 
to n = 0.02 given the sideband resolution, assuming the full finesse of 31,000, even 
for the 1.6 MHz mode. The silicon chip is mounted to holder also made from silicon. 
A compromise between high-Q m and mechanical stability is achieved by attaching the 
membrane chip only at three corners using Stycast 2850FT. The mechanical damping 
rate is measured in-situ by observing the mechanical ringdown lifetime with an optical 
probe at a wavelength where the cavity has low finesse. The realized values of Q m vary 
between 10 6 and 10 7 at 5 kelvin depending on chip mounting and membrane cleanliness. 
The matched thermal expansion of the chip and holder lowers the stress on the small 
epoxy joints. However, after thermal cycling from room temperature to cryogenic a 
few times, the epoxy joints are weakened and eventually the angular alignment of the 
membrane is compromised. 

The entire device is cooled using a 4 He flow cryostat from Advanced Research 
Systems Inc. Additional custom radiation shielding, including cold windows, is necessary 
to thermalize the membrane to the temperature of the copper cold finger (Fig. [2]). 
Alignment is monitored upon cooling the device by assessing the finesse of the cavity and 
the consistency of the optical mode position as a function of the membrane translation 
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Figure 2. Experiment layout. At the top is the laser and associated filter cavity. 
The three tunable acousto-optical modulators (AOMs) are for: The high-frequency 
branch of the cavity-laser lock, detuning the cooling laser frequency with respect to the 
locking laser frequency, and setting and maintaining the intensity of each beam. The 
locking laser is sent through an in-fiber electro-optical modulator (EOM) that applies 
frequency sidebands used for implementing a Pound-Drcvcr-Hall lock. Multiple beams 
are combined and sent directly via free-space into the cryogenic cavity. 



within the cavity standing wave. This analysis shows that we typically maintain sub- 
milliradian tilts between the components of the cavity. 

The cooling and probing of the membrane motion uses a 1064 nm Nd:YAG laser 
from Innolight (Fig. E^a)). While already a low-noise source, we additionally pass the 
light through a filter cavity (40 kHz linewidth) to remove intensity and frequency noise 
at our MHz frequencies of interest. Significant laser-cooling of our mechanical resonators 
requires careful attention to laser noise [12]. For a sense of scale, laser frequency noise 
of 1 Hz/a/Hz will result in a contribution in the output spectrum equivalent to a typical 
mechanical signal near the ground state. A single tone red-detuned from the cavity 
resonance is used for both cooling and detection via monitoring the transmitted intensity 
directly on a photodetector. We also make use of the orthogonally polarized mode for 
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locking the relative frequency of the laser and the cavity. At low frequencies we use the 
piezoelectric actuator attached to the mirror to stabilize the overall optical path length 
inside the cavity. High frequency noise is eliminated by servoing the laser frequency 
with an AOM. The overall bandwidth is 100 kHz, and aggressive filtering is used to 
eliminate response at the membrane mechanical resonance frequencies. Light from the 
two polarization modes is combined and split before and after the cavity on polarizing 
beam-splitters with cross-coupling of less than 10 -3 . 

With the membrane near the end of the cavity (Fig. [H(a)), the coupling between the 
membrane motion and the cavity resonance is more complicated than a membrane in the 
center of the cavity. The cavity resonance frequency, linewidth, and input/output port 
asymmetry all depend upon the placement of the membrane within the cavity standing 
wave, even without absorptive loss. These phenomena have been extensively numerically 
modeled in the work of Refs. [271 129] • A membrane placed near an end mirror can in 
fact have somewhat larger coupling than a membrane placed near the middle. However, 
operating at positions of enhanced coupling results in degraded finesse; hence, in our 
experiments we actually operate at the local coupling maxima where k is minimized. 

3. Results 

In this manuscript we present results from two different cavity devices. The first data 
shown in Fig. [3(a) uses a device optimized for cooling the (2,2) mode. As we increase 
the red-detuned input optical power, the mode is cooled and damped; the data show 
just the last amount of cooling. The mechanical quality factor is the relatively large 
value of Q m = 13.6 x 10 6 , and hence the mode is damped from an initial linewidth 
T m = 116 mHz to an optically damped T = 5.5 kHz (a factor of 47,000). 

If we expand the frequency range of the plot we see, however, that the overall 
relative transmitted intensity spectrum is not white (Fig. [3(b)). In these data there is a 
large peak to the left of the thermal peak of the (2,2) membrane mode. We attribute this 
peak to a thermally occupied mode of a part of the cavity mirrors/coatings or mounting 
structure (which we will refer to as "cavity mechanical modes"). Figure H]( a) shows a 
similar set of cooling curves for a second device focusing on (4,4) mode cooling. Here we 
see the cavity mechanical modes appear mainly as a weakly modulated spectrum at the 
highest optical powers (black curve). Similar features are evident near the (2,2) mode 
if a higher laser power is used to attain a lower shot noise level. We suspect the lower- 
ed modes in Fig. 11(a) arise from the fused silica mirror substrates based upon previous 
measurements in empty high-finesse cavities with similar mirrors [291130] , where at room 
temperature the modes were found to have a quality factor of 700. Particularly high-Q 
modes, such as the mode in Fig. [3] (Q m > 20, 000) we suspect may be from the single- 
crystal silicon substrate holding the membrane; our simulations of the free modes of the 
silicon substrate indicate there are indeed coupled modes around 1.6 MHz. In general, 
we have seen a variety of different cavity mechanical spectra depending on, for instance, 
substrate or mirror mounting techniques, coupling of the membrane motion compared 
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Figure 3. (a) Displacement spectrum inferred from cavity transmission near the 
(2,2) mechanical mode. Each curve is taken with a different cooling laser power, with 
cavity photon occupation ranging from 3 x 10 5 to 6 x 10 6 from top to bottom. For 
this device we operate with a cavity lincwidth of k/2it = 1.2 MHz and a detuning of 
A/27T = —1.6 MHz. (b) Relative transmitted intensity spectrum corresponding to the 
same data as in (a). The horizontal range has been widened to show the large peak 
on the left stemming from a mechanical mode of the cavity support structure. Note, 
detector noise that is a dominant contributor for the low power used in this dataset is 
not subtracted. Conclusions about phonon occupation for this device are displayed in 

Fig. m 



to the cavity end mirrors, or the membrane mechanical quality factor. 

The first step in the analysis of our data is a calibration of the membrane motion 
and an understanding of the membrane's bath temperature. We can extract the 
membrane coupling in three different ways: (1) Using thermally driven motion as a 
known displacement standard, (2) careful measurement of the cavity and membrane 
geometry including knowledge of membrane's position within the optical standing wave 
mode and a corresponding model of the expected coupling, and (3) the optical damping 
observed for a given intracavity photon number. In the Appendix we compare these 
models for the (2,2) mode data and find good agreement. The coupling for the (2,2) 
mode is found to be G/2ir = 1.9 x 10 16 Hz/m compared to an end-mirror coupling of 
u c /2ttL = 5.5 x 10 16 Hz/m. The (4,4) mode data presented here is not a particularly 
well-coupled device due to the final optical and mechanical mode overlap; for this device 
G/2tt = 3.9 x 10 15 Hz/m. 

To verify that the membrane is thermalizing to the expected cryostat temperatures, 
we look at Fig. E£b) that compiles data from the (4,4) device. We plot the 
integrated mechanical response from a Lorentzian fit to the data as a function of the 
mechanical damping (proportional to intracavity photon number) for three different 
cryostat temperatures. Here we plot only the larger response data where effects of 
cavity mechanical mode noise can be neglected. In Fig. H^c) we show the extracted 
bath temperature from these three datasets as a function of the measured cryostat 
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Figure 4. (a) Cavity cooling data for the (4,4) mode device. For this device we 
operate with a cavity linewidth of k/2it = 1.4 MHz and a detuning of A/2ir = —2.8 
MHz. As discussed in the text, based upon a calculation that includes the cavity 
noise, for these data we reach a minimum phonon occupation of n = 6; note a simple 
conversion between the intensity spectrum and motion (Eqn. [I]) yields n = 6 for 
the light blue curve (8th from the top). Detector noise is not a dominant source 
here, but the detection efficiency was low due to loss in the detection path. The 
small spike at 3.145 MHz in the lowest (black) dataset is electronic in origin, (b) 
Integrated mechanical response as a function of cooling power and cryostat temperature 
(points), along with fits to the data (lines). The three data sets correspond to cryostat 
thermometer readings of 4.9 (blue), 10 (green), and 15 (red) kelvin. The final point on 
the 4.9 kelvin data represents the light blue curve in (a), i.e. corresponds to n = 6 for 
a naive conversion, (c) Bath temperature extracted from the effective temperature in 
(b) as a function of cryostat thermometer temperature. The vertical error bars are the 
uncertainty extracted from the uncertainty in G and mass, and the horizontal error 
bars represent the systematic uncertainty of the thermometer. 



temperature (taking into account slight variations in mechanical and cavity parameters 
with the changing bath temperature). We find a linear trend that extrapolates to zero 
within the uncertainty, indicating that the membrane mechanical mode does indeed 
equilibrate to the cryostat cold finger temperature. 

To start we can use our extracted coupling G to apply a naive conversion between 
the transmitted intensity spectrum and the mechanical motion of the mode of interest. 
When cavity mechanical motion and radiation pressure shot noise can be neglected, 
we can convert the relative transmitted intensity spectrum Sj(u) to a (one-sided) 
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displacement spectrum S z (cu) using the following relation: 

Sz{u) = ~7^~ |II(a;)| 2 G 2 ' (1) 

Here U(co) = Xc(u) — u ) where Xc(u) = K /2-l(u+h) * s ^ ne res P onse function for a 
cavity with linewidth n (energy decay rate) and detuning of the laser from the effective 
cavity resonance A. The extracted S z (u) yields a thermally driven, optomechanically 
damped Lorentzian on top of a floor coming from optical shot and photodetector dark 
noise (Fig. |3]^a) for example). We extract the effective temperature of the mode by 
fitting the data to a Lorentzian. We use the area under the fit to determine the mean 
square displacement, which is proportional to the effective temperature, and the effective 
phonon occupation n. Applying this procedure to the (2,2) mode data in Fig. [3]^a) we 
find the apparent effective phonon occupation displayed as the red squares in Fig. |5^d), 
which would indicate the thermal component of the mechanical motion is brought to 
1.4 quanta. 

As discussed above, an important limit to laser cooling is classical noise in the 
relative frequency of the cooling laser and cavity resonance frequency. Such effects 
include laser phase noise [311 [32j [33] and mechanical and thermal noise in the cavity 
mechanical modes [3U El [30]. This relative frequency noise is converted to intensity 
noise in the cavity for off resonant laser drives, such as cooling tones. The intensity 
noise applies a radiation pressure force to the mechanical oscillator leading to an extra 
displacement and greater effective temperature. Further, extracting the temperature 
from the transmitted intensity spectrum in the presence of classical noise becomes 
complex [HJ[321[I2]. 

As a first study of the effect of the classical noise on measurements of n, we present 
a more thorough analysis of the data of Fig. [3]^b). Here, we model the cavity mechanical 
mode noise as dominated by a single thermally-occupied mode. We then determine the 
corresponding effect on the (2,2) mode. Importantly, we assume that the observed cavity 
mechanical motion does not directly drive the membrane mode, but rather the only 
coupling is through the intracavity field. If this assumption were not true, it is possible 
the entire output spectrum could be attributed to membrane motion. Further, we 
assume the effective mass of the cavity mode is much larger than that of the membrane 
mode. The physical mass of expected cavity modes is greater than 10 6 times larger than 
that of membrane modes. Hence, even without complete knowledge of the coupling to 
the cavity mode, this is an excellent assumption. Thus, the cavity mechanical mode 
motion, to a good approximation, is unaffected by the optomechanical interaction. 
In the Appendix we present a calculation for this situation. Our model allows us to 
estimate the displacement spectrum of the membrane mode from the measured cavity 
transmission spectrum, when the simple relation between transmission and displacement 
spectra (Eqn. [1]) no longer holds. 

We estimate the functional form of the noise spectrum as a Lorentzian centered at 
the position of the large cavity mechanical mode peak. Because the noise is concentrated 
away from the membrane mode resonance, we may assume (at least for the weakly 
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Figure 5. (a) Si(w) designed to model the observed relative intensity spectrum for 
the data in Fig. |3jb) at an intracavity cavity photon number of 3.3 x 10 6 (teal curve 
7th from the top in Fig.[3jb)). (b) Expected contributions to the relative transmission 
spectrum due to the cavity mechanical modes, <Ss(w) for the same parameters as in 
(a), (c) Modeled displacement spectrum of the (2,2) mechanical mode (solid green). 
The dashed red line shows the displacement spectrum expected if we were to ignore 
Si(uj). (d) Extracted integrated motion (converted to a phonon occupation) for the 
(2,2) membrane mechanical mode as a function of cooling laser power. Red squares 
show only the contribution of the thermally driven Lorcntzian peak, which correctly 
predicts a small on-resonance motional spectral density, but gives an erroneously low 
phonon occupation. The green diamonds include contribution from red squares and the 
motion induced by the cavity mechanical mode motion. The curve in (c) corresponds 
to the point with a 3 kHz mechanical damping rate. 



damped data) that the area in the transmission spectrum due to the cavity mechanical 
mode is largely independent of the membrane motion, and is equal to Si(u), the 
transmission spectrum without the membrane mode of interest present. In the Appendix 
it is shown that this extra intensity noise causes a radiation pressure force that drives 
the membrane mode to an additional displacement given by: 

S *M W ) - p V (2) 



In this expression, g = GZ zp where Z zp = w 2m ^ is the zero point motion for an 

oscillator with frequency u m and effective mass m, and iV is the intracavity photon 
number. The function J\f(u) represents the optomechanically modified mechanical 
response and is displayed in the Appendix. 
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Curves displaying an example of the modeled spectrum Sj(u) and the cavity 
mechanical noise contribution Si(u) are shown in Fig. [5(a) &(b). The inferred value 
of S z (oj) for these parameters is plotted in Fig. [5(c) (green line). Since the majority of 
the noise peak in Sj(uj) comes directly from Si(u), the height of the cavity mechanical 
mode peak relative to the thermally driven motion in the displacement spectrum is much 
smaller than the corresponding ratio in the transmission spectrum. The total estimated 
displacement of the membrane mode is shown in Fig. [5](d) as the green diamonds. 
Specifically, we plot the integrated motion (z 2 ) for curves like that shown in Fig. [5(c), 
converted to an equivalent phonon occupation. As the laser power increases, the 
contribution of the cavity mechanical mode noise drive increases because the intensity 
fluctuations grow in proportion to the optical intensity. Thus, a minimum phonon 
occupation of n ~ 5 occurs at an intermediate intracavity power. And as noted 
earlier, without this cavity mechanical mode (or if this mode were much farther away in 
frequency space) we would achieve n approaching unity, as displayed by the red squares. 

Next we consider the temperature achieved for the (4,4) mode data of Fig. 11(a) 
taking into account the cavity noise. Here our best estimate is obtained by modeling 
the cavity mechanical noise spectrum as a white floor. For these data, shot noise 
for the highest-power set (black curve) corresponds to a relative intensity noise of 
6 x 10~ 15 and hence we would place the classical cavity noise floor at ~ 4 x 10~ 15 . 
First, it is useful to consider the magnitude of motion this value corresponds to; as 
discussed above, we suspect in this case the cavity mechanical motion is dominated by 
the mirror fused silica substrates. Using the G values discussed above, Si/ I 2 = 4 x 1CT 15 
corresponds to membrane motion of 0.9 x 10 -17 m/yTlz or a real end-mirror motion of 
0.7 x 10~ 18 m/vTlz- Assuming a scaling of the mirror substrate motion with A/^bath, this 
cavity mirror motion is within a range expected for thermally occupied modes of fused 
silica substrates based upon previous measurements of room temperature cavities [29J . 
Finally, we apply the same model discussed above and in the Appendix to these data. 
Again we assume that there is no physical coupling between the cavity mechanical modes 
and the membrane; since we assign the cavity mechanical modes in this situation to the 
mirror substrates, we believe this is a very good assumption. We use the estimated flat 
spectrum at 4 x 10~ 15 as Si{uj)/I 2 and the data of Fig. Ufa), and we find the phonon 
occupation reaches a minimum value of n = 6 as a function of cooling laser intensity. 

In conclusion, we have presented two experimental data sets in which we have cooled 
SisN4 membrane modes, as a conservative estimate, within a factor of 10 of the quantum 
mechanical ground state. The dominant uncertainty in our temperature measurement 
comes from complex classical noise spectra added at the few phonon level by other 
mechanical modes within the cavity. We expect with minimal additional optimization 
of our cavities we can remove deleterious effects of these additional cavity modes. 
Already we have demonstrated the kind of system parameters required for ground state 
cooling; for example a fruitful set of parameters would combine the coupling and quality 
factors achieved for the device in Fig. [3] with the absence of the isolated high-Q peak 
in Fig. HI We believe engineering the silicon substrate to advantageously define the 
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relative frequencies of the membrane and substrate modes will allow us to consistently 
achieve desired parameters in future designs [35]. It has also been shown that relative 
frequency noise of the cooling laser and cavity resonance frequency can be removed 
via active feedback schemes involving higher-order cavity modes [30]. Further, final 
studies of these devices will likely be conducted at even colder cryogenic temperatures 
compatible with superconducting circuits; here thermal occupation of cavity mirrors 
and the substrate should be reduced. 

Importantly, thus far we have not observed physical absorption heating of our 
devices that leads to a significant increase in Tb at h- This fact, combined with the efficient 
detection that these devices can afford, makes this system promising for the observation 
of the radiation pressure shot noise. Even with the presence of cavity frequency noise, 
shot-noise from an sufficiently-strong optical tone placed on the cavity resonance should 
be efficiently transduced to the membrane motion. 
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Appendix A. Calculation of thermo-mechanical noise in an optomechanical 
system 

To analyze our optomechanical system we consider a standard optomechanical 
Hamiltonian with the addition of cavity mechanical modes of the cavity structure 
indexed by i: 

H = hjjj m c^c + hw c a)a + HGZ zp (c + c*)a)a (A.l) 
+hJ2 (^i&i&i + GiZ^ibi + b\)a!a\ +H K + H T 

i 

Here c is the annihilation operator for the mechanical mode of interest (the membrane 
mode), with oscillation frequency u m and harmonic oscillator length Z zp . a is the 
annihilation operator for the cavity mode, at frequency u c . G is the optomechanical 
coupling constant, and we define a single photon coupling rate go = GZ zp . bi is the 
annihilation operator for the ith cavity mechanical mode with frequency Ui, harmonic 
oscillator length Z zpi , and optomechanical coupling Gj. Additionally, the term H K 
represents the input and output optical coupling of the cavity with total cavity decay 
rate k = k l + n R + km- The contributing decay rates stem from the input port («&), 
the output port (kr), and the internal loss (/tint). The term represents the thermal 
drive on all of the mechanical modes. 

The optomechanical interaction can be linearized and fast oscillations at the optical 
frequency accounted for by defining a(t) = (a + d{t))e~ luJLt . a = 'fj^-i'L is the large 
classical amplitude of the intracavity field and d(t) represents small fluctuations about 
this value. d m is the coherent state amplitude of the input laser field driven at frequency 

U) L . 

From the Hamiltonian in Eqn. ( lA.lj) . we derive a set of Heisenberg-Langevin 



equations of motion and transform them into the frequency domain using the Fourier 
transformation f(u) = J^e 1 ^ f(t)dt, p(u) = j^e^ 1 f\t)dt. 



d(u) = Xc(w) (-iaGz(u) - ig^GiZiju)) + ^(w) + vW6ntM + v^fflM 

zp 

z and Zi represent the small displacements of the mechanical modes about their 
optomechanically shifted equilibrium positions z and Zi, such that Z zp (c + c^) = z + z 

and Z zP:i (bi + b\ \ = z t + Zi. We use here the mechanical susceptibilities Xm(w) = 

r m /2-l(u,-u m ) and = Ti/2 J (U - Wi y Xc{u) = k/2 J {a+uj) is the cavity susceptibility, 

where A = u c — — Gz — ^GiZi is the detuning of the laser input frequency from 
optomechanically shifted cavity resonance frequency. The operators £_Le~ iaJi *, ^Re~ luJLt , 
and ^mtC" ;iWL * are Langevin noise operators representing vacuum fluctuations entering 
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the cavity from the input, loss, and output ports (see for example [SU EZ])- V an d T)i 
are the Langevin noise operators representing the thermal and vacuum noise driving the 
mechanical modes. To simplify the equations of motion we drop small terms of order 
d 2 , dz, and dzi. Additionally, we work in the limit where G 2 Z 2 pi a*d/ k <C n th i Yi where 
n t h,i = fc ^ ath and T bath is the bath temperature of the modes. In this limit the optical 
drive does not perturb the mechanical state of the cavity mechanical modes. However, 
signatures of the cavity mechanical modes are still imprinted on the optical mode. 

We can then solve the equations of motion and calculate the mechanical 



displacement spectrum S. 
(z(-u)z(u)) 



(2), 



Af(-u)Af(u) 



r 



-u)z(u)}. 



+ 



n t h 



Xm(w)| 2 



Xm(-W)| 2 
-C)| 2 



(A.2) 



+ 4^ 2 (a*a) 2 |n(u,)| 2 (6f(-u)6f(u)) 



where Af(u) 



— i2u m g 2 a* aXl{uj) and n(oj) 



make use of the operator expectation values 



,Cint(-^)Cint( 



We 



1, and {rj{-uj)r] ] '(u)) = n th + 1, 



n t h- 



The 



first term of Eqn. ( 1A.2I) contains the thermal motion of the membrane. The second 
term includes the membrane motion induced by radiation pressure from optical shot 
noise. The effect of the cavity mechanical modes is seen in the third term, through 
(5f(—u))5f(u))) = J2G 2 (zi(—u)zi(u)), which represents the noise spectrum of the cavity 
frequency shifts induced by the cavity mechanical modes. 

We next compute the spectrum of intensity fluctuations for light directly detected on 
a photodetector at the output port of the cavity. Let Sjf = ((/(— w) — 7) — I)) 

be the two-sided power spectrum of the detected photocurrent where I(t) = 
efuv L Ual nt (t)a out (t) + (1 - e)a\ l (t)a n (t) + I d (t), with mean value I = (I). TZ = 
is the photodetector sensitivity where q e is the electron charge, e is the detection 
efficiency, a n (t) = £, n {t)e~ luJLt is a Langevin noise operator representing vacuum 
fluctuations entering the detector through the loss port associated with the detector 
inefficiency, and I d is the photodetector dark current. The output optical field a out (t) = 
(flout + d out (t))e~ luJLt is evaluated via the input-output relations: a out = VK^a and 



R 



Had. 



s 



(2) 
// 1 



p 



K R [a*a) 



— * 1 
a { 



K R d(- 



X 



1 



/ ~KR~d(LJ 
1 



- £r{-u))) + a(^KR~d ] (-u 
(I d (-u)I d (u)) 



-a,)) 



e K R a*a 



I 2 

c 



e K R ara 



+ 



[I d (-u)I d (u)) 
P 
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10J 



y/ K R 

^ m = -iGIL(u)z(u) 

Vi = -iU(u)5f(u). 
We can thus compute the functions 



[9 q {-u>)* q { U )) 



1 



i2 1 



(V m {-u)V m {u;)) = G 2 \U(u)\ 2 (z(-u)z(u)) 

(%(-u)%(u)) = \U(u)\ 2 (5f(-uj)5f(uj)) 
and the nonzero cross terms are 

($ q (-u>)y m (uj)} + (& m (-u>)y q (u)} = -Au m g 2 lm 



U(u) 



Xc 



-Auj m a* agl\m. 



-LU 



|n(u;)| 2 (5/(-u;)5/(u;)). 



The photodetector signals we record in the experiment are one-sided power spectra 
Sj(uj) = S$(-w) + Sffitu). Similarly we define S l (u)/P = 2 <*i(-w)*i(w)), the 
one-sided spectrum of transmitted intensity noise due to the fluctuations of the cavity 
mechanical modes, and from this we infer S z (u) = (z(—cu)z(lo)) + (z(uj)z(—u)) the 
one-sided displacement spectrum. 

If the contribution of the cavity mechanical modes is small (Si(u) — > 0) and the 
radiation pressure shot noise is small, then the displacement spectrum can be easily 
inferred from the transmission spectrum. 

1 SAu) 



S z ,8f^o( UJ ) + ^'noisefloor(k') 

Where S no i SC f[ oor (uj) = G 2^^ 



G 2 |IT(u;)| 2 P 

2 K Ja*a + (Id(—u)Id(u)) /I 2 ^ is the detection noise floor 
due to optical shot noise and detector noise. If frequency noise is not negligible 
and Si(u) can be measured empirically then the contribution to S z (u) from the 
optomechanically transduced fluctuations of the cavity mechanical modes can be 
computed from Eqn. ( 1A.2I) . 



S Z)5f (uj) _ Au 2 m gl{a*a) 2 Si{u) 



(A.3) 



\M(uj)\ 2 P 

Note a* a is the intracavity photon number N. From S z (u) the root mean square 
displacement of the membrane mode can be computed. 



z 2 p 



S z (u) doj 



2U+ 5 
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We can compare our final expression for the impact of mechanical modes within 
the cavity (Eqn. IA.3j) to the result expected for equivalent frequency noise on a laser at 
the input port of the cavity. We find our expression can be translated into a formula 
equivalent to the laser frequency-noise result derived in, for example, Ref. |32j . 

Appendix B. Calibration of membrane motion 

Here we present the three different methods we use to calibrate the mechanical motion 
and corresponding temperature. We compare the methods by assessing the membrane 
coupling extracted in each case for the (2,2) data presented in Fig. [3j based on these 
analyses we deduce an uncertainty in G of 5%. 

Appendix B.l. Bath temperature and optical damping 

The thermal motion of the mechanical resonator is equated to the effective temperature 
T via (z 2 ) = ^jr- Optomechanical cooling theory in the large damping limit tells 
us T = Tbath--j?S where T can simply be determined by the measured linewidth of the 
optically damped mechanical resonator. We compare (z 2 ) as determined via integrating 
the spectrum J °° S z (u)^ to that from fc ^| th ^. T bath = 4.9 K is input based upon the 
cryostat thermometer, where the accuracy of this measurement was motivated via the 
trend in Fig. |U[c). S z (u) can be extracted from the measured intensity spectrum via 
Eqn. [T]with one free parameter G. 

The inputs to this calculation are: k/2tt = 1.2 MHz, which is determined by a 
measurement of the cavity linewidth in ringdown for the position of the membrane 
during the measurement. A/2n = 1.6 MHz, which is determined from the detuning 
set with respect to the on- resonant locking light, accounting for a birefringent splitting 
of 0.4 MHz. The mechanical frequency u m /2n = 1.575 MHz. The mechanical quality 
factor Q m = 13.6 x 10 6 . The effective mass m = pd 2 t/4, which is determined using 
p = 2700 kg/m 3 [38]. (However, given the range of values in the literature for LPCVD 
Si 3 N 4 [29] we would place a 10% systematic uncertainty on the mass. Note, the mass 
uncertainty is relevant for G but not go [39].) 

This comparison tells us G/2tt = 1.8 x 10 16 Hz/m. 

Appendix B.2. Membrane- at-the- end model 

With knowledge of the position of the membrane within the cavity standing wave we can 
model the expected coupling. In each iteration of the experiment, we scan the membrane 
within the cavity to sit at the position for which the cavity linewidth k is minimized [29J . 
With knowledge of the cavity and membrane parameters we can calculate the expected 
coupling at this operating position we refer to as Z min . Z min is registered within the 
standing wave ~ 0.9 mm from the flat mirror of the cavity. The cavity has an overall 
length of 5.1 mm, which is measured via the ratio of the transverse mode spacing to the 
free spectral range of 29.4 GHz. Based on direct ellipsometry measurements of our film 
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we base our calculation on a t = 40 nm thick membrane with index n = 2.0. Hence, at 
Z m in we predict du c /dz = 2tx x 2.9 x 10 16 Hz/m. 

We then apply a correction based upon the measured transverse mode overlap 
between the mechanical mode of interest and the TEM o cavity mode (an example of 
such a measurement is shown in Fig. [1]). To account for the mode overlap we calculate 
G = rj^ using 00] 



JJ dxdyl(x,y)a mn (x,y)/a 



(B.1) 



where a mn = ao sm(mirx/d) sin(n7ry/d) and I(x,y) is a normalized intensity function 
I(x,y) = ?rw 2 w exp(— 2(x — xo) 2 /w^)exp(— 2(y — yo) 2 /Wy). We measure the location of 
the (2,2) mode while cold to be at the coordinate position (x ,y ) = (108,99) /im. The 
size of the mode in x,y is measured to be w^ = 92 /xm and w y = 88 /xm. Hence 77 = 0.67, 
and G/2n = 2.0 x 10 16 Hz/m. 

Appendix B.3. Optical damping and intracavity photon number 

The optical damping V we observe for a given intracavity photon number is also a 
measure of G. To infer the intracavity photon number N from the measured output 
flux of the cavity we understand the asymmetry of the cavity and the internal loss. 
Usually we orient the cavity according to Fig. [1] where the membrane is at the output 
(right) side of the cavity, providing the most signal at the output. However, for these 
particular measurements we happened to measure out the less-transmissive port of the 
cavity, i.e. the membrane and flat mirror were actually at the input (left) side of the 
cavity. For a measurement of the photocurrent I at the output port the photon number 
is given by 

— Ptw\i K> 1X1 . . 

]V= ou L _ = (B.2) 

Fuj c k k r k r q e e 

for a low-reflectivity membrane [29], where e = e^e^ is a combination of two efficiency 
factors. In our setup for these measurements, the detector efficiency is e& = 0.87 and 
the propagation losses from the cavity output to the detector are given by e p = 0.88. 

To find we must understand all the contributions to k = kr + k,l + K- mt . K- mt 
is dominated by clipping of the transverse mode on the silicon frame that results from 
imperfect alignment when cooling to cryogenic temperatures, and hence for modeling G 
and k we consider this a loss that is independent of membrane position. At the position 
that corresponds to the minimum value of k (Z min ), the theoretical analysis of the three- 
element cavity described above indicates /t m ; n = 0.79 MHz. With a well-aligned cavity 
we often achieve this value; in this particular cavity at room temperature we achieve 
ftmin = 0.85 MHz, and when cooled down we find K m i n = 1.17 MHz. This indicates an 
internal loss contribution of K- int = 0.33k. The ratio between kr and kl can be calculated 
from the asymmetry of the cavity, which can be determined from a calculation of the 
resonant reflection R and transmission T at our operating position of Z m [ n . We use the 
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expression kl/kr = — "r~ ^° ^ nc ^ k l/kr. = 1-9 [29]. Thus /€r = 0.23k taking into 
account n int . 

We then compare the measured optically damped linewidth T to the expected 
calculated T given by [HJ W2\ 



This gives us G/2tt = 1.9 x 10 16 Hz/m. 
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